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Abstract: In this article, a class of Fourier Integral Operators which converge to the unitary 
group of the Schrodinger equation in semiclassical limit e — > is constructed. The convergence is in 
the uniform operator norm and allows for a error bound Cjve^"'"'^ for any integer A'^ and extends to 
Ehrenfest timescaleswith bound Cjve^'*"^"'' where p can be made arbitrary small. In the chemical 
literature those approximations are known as simple Initial Value Representations. 

1 Introduction 

We study approximate solutions of the semiclassical time-dependent Schrodinger equation 
iej^r{t) = -'-^^i''{t) + V{x)rit), r{0)^i'l^L\W\C) (1) 

2 

in the semiclassical limit e — > 0. The operator H' :— — -f- V{x) on the right-hand 
side of ([T]) is the so-called Hamiltonian, a self-adjoint operator on _L^(M'^) (under suitable 
assumptions on the potential 1^). It is well-known that the solution of ([T]) can be written as 

where the group of unitary operators e^i^ * is defined by the spectral theorem. 

The semiclassical parameter e may be thought of as the quantum of action h, but there 
are also situations, where e has a different meaning. One example is provided by Born- 
Oppenheimer molecular dynamics, where equation ([T]) describes the semiclassical motion of 
the nuclei of a molecule in the case of well-separated electronic energy surfaces and e is the 
square root of the ratio of the electronic mass and the average nuclear mass. In this case, the 
e in front of the time-derivative in ([T]) is due to a rescaling of time t = t/e. This particular 
choice, the so-called "distinguished limit" (see [S]), produces the most interesting results in 
the semiclassical limit e ^ 0. 

To formulate our main result, we introduce the following class of Fourier Integral Opera- 
tors (FIOs): 

I^(/t*;w)</)(x):=-i-^ / ei-'^^^--'y'^^uix,y,r,)q>{y)d7jdy, (2) 

where 

• K*{q^p) ~ (g,p),S'' iq,p)j is a C^-family of canonical transformations of the 
classical phase space T*R'' = R'' x R'^, 

• S'^ {QtP) is the associated classical action (see Definition [3] for a precise definition) 
which reduces to 

t 



{q,p) = 



dT 





when K* is the Hamiltonian flow associated with a Hamiltonian function h, 
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• the complex- valued phase function is given by 

• and the symbol u is a smooth complex-valued function which is bounded with all its 
derivatives. 

The central results of this paper reads 

Theorem A For fixed t, under the assumption that the canonical transformation At* is 
hiLipschitzian with hounded derivatives, the operator I'^ [k^ ; u) is continuous S{M.'^) — > iS(K'^) 
and extends to a continuous operator L^(R'^) L'^{M.'^) with 

Theorem B Let e~e^^* be the propagator defined by the time- dependent Schrddinger 
equation ^ on the time-interval [— T, T] with subquadratic potential V E C°°(R'',R), i.e. 
sup^gjjd \d"V{x)\ < oo for all a e N'^ with \a\ ^ 2. Then 



sup 

te[-T,T] 



^ C{T)e, 



where k* = (X" , ) and u are uniquely given as 

• the flow associated with the classical Hamiltonian h{x,£_) = ^I^P + V{x) 



-VV (x^'ig^p) 



X-^{q,p) = q 



the solution of the Cauchy-problem 



dt 



uit,y,ri) 



-tr 



1 

-1 

2 

1. 



u{t,y,r]) 



M(0,y,?7) 
The C'^^'^ -valued function 

yo (y , 7,)) = S,f (y , 7?) - ^xf {y,r,), 
depends on elements of the transposed Jacobian 



of Kf with respect to {y,ri). 



The equation for u is easily solved. Its solution 

u{t,y,7]) = (dcty(F''\y,r^) 



(where the branch of the square root is chosen by continuity in time starting from t ~ 0) 
mirrors the so-called Herman-Kluk prefactor as presented in [17] . The separation of the 
results in two steps is made on purpose to emphasize that the first L^-boundcdness result 
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is a key argument into the proof of the second result. By the way, Theorem B is just a 
simphfied version of our main result in Section [5) More precisely, Theorem |3 will essentially 
add two central aspects. First, for the Ehrenfest-timescale T(e) = Ct log(e~^) the result 
still holds with a slightly weaker bound. Second, the error estimate can be improved to , 
where N is arbitrary large by adding a correction of the form J2n=i ^"''^n to u. As u, the w„ 
are solutions of explicitly solvable Cauchy-problems. Finally the opportunity to extend this 
result to general pseudodiffcrcntial operators will be discussed in Remark [5j 

Whereas there is an abundant number of works on Fourier Integral Operators in the 
mathematical literature, only few of them discuss the relation between FIOs and the time- 
dependent Schrodinger-equation. The first works which apply FIOs with real-valued phase 
function to this problem are [H] and [TT]. In this case one has to deal with the boundary 
value problem 

Given a;, y G K'', find p such that = x. 

To get uniqueness for its solution one has either to restrict to short times t or to impose 
very strong restrictions on the potential. The same problems are met in [7], where Fujiwara 
applies a related class of operators without integral in the oscillatory kernel to the Schrodinger 
equation to justify the time-slicing approach for Feynman's path integrals. 

The avoidance of this problem is the major advantage of complex- valued phase functions. 
In the non-semiclassical setting, Tataru shows in [18] that the unitary group of time evolution 
is an FIO with complex-valued phase function (different from ([3]) because his kernel consists 
of an integral over the whole phase-space space in contrast to the momentum integral we 
proposed). He also establishes that the use of constant symbols leads to a parametrix for the 
non-semiclassical Schrodinger equation. In the semiclassical setting, [5], [13] and [IT] treat 
FIOs of the same type (phase-space integral kernel) with closer connection to Gaussian wave 
packets known in the chemical literature as Initial Value Representations. 

The class of operators defined in ^ is used in the works [T31 and [3] for the construction 
of approximate solutions of the semiclassical time-dependent Schrodinger equation. However, 
these works only allow compactly supported symbols, which enforces the truncation of the 
Hamiltonian at least in momentum. In particular, in [13j the comments of Theorem 2.1 
mentioned that the dependence of a constant analog to our C{T) on 

• the size of the compact support, 

• the length of the time interval T and 

• growth of a possible subprincipal symbol for iJ^ 

was not controlled. The result presented here proposes a (partial) answer to those questions. 

In contrast to the mathematical literature connecting time-dependent Schrodinger equa- 
tion and Fourier Integral Operators, there is an abundant number of papers in chemical 
journals on this topic. The reader interested in those aspects can consult |10j (and references 
therein) which gives a good review of the different approximations used by chemists. 

Organisation of the paper and notation 

The paper is organised in the following way. Section [5] will introduce the class of FIOs 
with complex quadratic phase and constant spreading matrix O. In Section [3l we will prove 
their continuity on the Schwartz space 5(M'^) and explain how to push the dependence of 
the symbol on the variable x to higher order in e. Section U] will adress the problem of the 
L^-boundedness for those FIOs. The extension of all those results to non-constant spreading 
matrix Q{y,T]) will be hinted in Section [5l Finally, all these results will lead to our main 
result, which we will state in Theorem |4] of Section [6l 

We close this introduction by a short discussion of the notation. Throughout this paper, 
we will use standard multiindex notation. Vectors will always be considered as column 
vectors. The inner product of two vectors a, £ R'' will be denoted as a • & = '^j^j ^^'^ 
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extended to vectors a,b & by the same formula. The transpose of a matrix A will be A'' , 
whereas A* := A^ denotes the adjoint and finally ej will stand for the jth canonical basis 
vector of or C''. 

For a differentiable mapping F E C^(R'^,C'^), we will use both {dxF){x) and Fx{x) for 
the transpose of its Jacobian at .t, i.e. {{dxF){x))jk = {Fx{x))jk = {dxjFk){x). This leads 
to the identity dx{F ■ G) = GxF + F^G for F, G G G^{R'^,C'^). The Hessian matrix of a 
mapping F € C^(M'',C) wiU be denoted by Hess2;i^(x). 

For the sake of better readability of the formulae, we will be somewhat sloppy with respect 
to the distinction between functions and their values. As a crucial example, we will write 
{x — X'^{y, r]))v for the function (x, y, ?/) i-^ (x — X'^{y, r]))v{x, y, rf). 

Finally, we remind the normalization for the Fourier transform we will be using 

^ £ e-^>(.)d. and = £ e^^^m^i. 

2 Semiclassical Fourier Integral Operators 

Our FIOs involve two fundamental objects: their symbol u and a canonical transformation 
K. We will first give more precision to the restrictions wc put on those objects. 

Definition 1 (Symbol class) Let d = {dj)i^j^j € N"^ , u[z\, . . . , zj) a smooth function of 
C°°(R''i X • • • X M"^'; C^) and m = (mj)i^j^j € R-' . We say that u is a symbol of class 
S'[m; d] if the following quantities are finite for any k ^ 

M™[u] := max sup 



where (z) := ^/l + \z\^ . 

We extend this definition to any rrij G M := {— oo} U R U {+00} by setting, for instance 
with non-finite mi, 

5[(+oo,TO2, ■ • ■ ,'Tij);d] = y ^[(mi, . . . ,mj);d] 

mieR 

and 

5[(-cx),TO2, . . . ,mj);d] = Q 5[(mi, . . . , mj); d] 

and so on. 

To fix notations, we also recall the definition of a canonical transformation and the link 
with symplcctic matrices. 

Definition 2 (Canonical transformation) 

Let K{q,p) be a smooth diffeomorphism from R'^ x R'' into itself decomposed into posi- 
tion/momentum variables through K{q,p) = {X'^{q,p),'E,'^{q,p)). We represent its differential 
by the following Jacobian matrix 

^ - 1, s^(g,p)t E;{q,py ) ■ 

K is said to be a canonical transformation ifF'^{q,p) is symplcctic for any {q,p) mR'^xM'^ 
i.e. 

[F'^{q,p)]^JF''{q,p) = J where J := 
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As discussed in |15] and [17] , one can associate to a canonical transformation a real- valued 
function that plays the role of the Lagrangian action integral for Hamiltonian flows. 

Definition 3 (Action) Let K{q,p) ~ (X'^{q,p),'E.'^{q,p)) be a canonical transformation of 
R'' X M''. A real-valued function S'^ is called an action associated with k if it fulfills 

S^iq,p) = -p + X^iq,p)E^{q,p), S;{q,p) = X^{q,p)E''iq,p). (5) 

Remark 1 The function S"" always exists and is uniquely defined up to an additive constant. 

We specialize here the notion of diffeomorphism of class B as presented by Fujiwara in [7]. 

Definition 4 A canonical transformation k of x M.'^ is said to be of class B if F'^ G 

^[O; 2d] and for any k^O, we set MllF'']. 

Remark 2 The subset of canonical transformations of class B is a subgroup (for composi- 
tion) of bilipschitzian diffeomorphisms ofW'' x R'' i.e. if k is a canonical transformation of 
class B, then there exist two strictly positive constants and such that for any {qi,pi) 
and {q2,P2) in M'' x R"^ 

CK||(g2,P2) - {qi,pi)\\ < ||k(92,P2) - K{qi,pi)\\ ^ C^||(g2,P2) - (gi,pi)||- (6) 

From now on, all canonical transformations considered are assumed to be of class B and 
e will denote a small parameter such that < £ ^ 1. 

Wc turn now to the definition of the main object of this article: semiclassical FIOs with 
quadratic complex phase. 

Definition 5 (FIO) Let O be a complex symmetric matrix with positive definite real part. 
For u E 5[+oo; 3d] and ip E S{M.'^; C), we define the action on (p of the semiclassical FIO 
associated with k with symbol u as the oscillatory integral (see Avvendix V^) 

[I'{K-u-Q)ip\{x) :=7;r^ / eJ*^(--^''''^)Kx,y,,7Mj/)]dr;dy 

where is a complex-valued phase function given by 

^^[x, y, 77; 9) = S'^iy, 77) + 77) • (x - X-{y, r,)) + '^{x - X^{y, r,)) ■ e(.T - X^iy, 77)). 

Remarks 3 

1. For any u G 5'[-|-oo; 3d], it results from the usual machinery of oscillatory integrals that 
the operator [k; u; Q) is continuous from S{W^) into its dualS'{W^). 

2. The normalization is adjusted so that we have X^(Id; 1; O) — Id. 

3. With the rescalings 

• T^ip{y) £'^l^'p){yfey) (which is unitary on i^(M'')j, 

• K^'^\y,if) := K,{^y ^ ^r\) I yfe (which preserves the symplectic structure), 

• u'-''^{x,y,il) := u{y/£x,y/ey,y/eT]), 

we have 5"=*"' (7/, 77) = S'^ {^y , ^/erj) / e , 

(x, y, 77; 6) = ^''(Via;, Viy, ^/i77; Q)le 

and 

i^K-u-Q) = (r|)*xi(^(^^7.(-);e)T|. (7) 

4. One might connect those FIOs to the one defined in |15j and |17j by taking the limit 

+00 in the latest with an appropriate renormalization. However, it turns out that 
this limit process does not allow to transfer their more interesting properties ( essentially 
L"^ -boundedness) to X"^ {k; u; Q) . 



5 



3 S Continuity and Integration by Parts 

The continuity obtained for the FIOs of [iS] and [l7] on the Schwartz space 5(M'') can easily 
be obtained for T^{k; u; 0) by the same type of argument. 

Theorem 1 If u G S[+oo;3d], then I'^ {k; u; Q) sends S{M.'^) into itself and is continuous. 
Moreover, for finite rrf , and m^, the map 

u G S[{m'' ^m"); {d,d,d)] ^ (X"(K;u;e) : S{m.'^) 5(M'^)) 

is continuous. 

The proof follows from various integration by parts. We do not give a detailed proof here 
but remind the main tools. 

We compute the derivatives of with respect to the variables x, y and rj 

iO \ / E-{y,r,) \ 

S^-*X«e] iy,v) \\ V ■ (8) 

S--zX-e] (y,7?) / V x-X-iy,r^) J 

In order to integrate by parts, we need the invertibility of this 3d x 3d matrix so we begin 
by establishing invertibility properties of 

y{F; 6) := ( I ) ( ^'^^ ^ ^ ~ iB^Q (9) 





for a matrix F with block decomposition 



A B 
C D 



Lemma 1 If F is symplectic, y{F;Q) is invertible and if k is of class B, y{F'^{-, ■);&) ^ 
is in the class S[0;2d\. 

Proof A straightforward computation shows that 

y{F; 6) {my'^ y{F; ey = v{F; e)*v(F; e) 

where 

Hence the invertibility of y{F; O) as the kernel of V{F; 8) is reduced to {0}. 

If K is of class B, then y{F'^; 8) is clearly in S[Q; 2d] so the formula of the inverse with 
minors shows that it is enough to prove a bound from below for the determinant. y{F; 8) 
is linear continuous in F and, as k is of class B, -F" has relatively compact image, hence 
I det3^(_F; 8)| reaching its minimum on this compact, minimum which can not be zero from 
above. □ 

Thus, integrating by parts with respect to rj allows to convert any polynomial growth in 
X to the variables (y,f?) through the decomposition x = X'^{y,rf) + (x — X^{y,ri)). More 
precisely, the removal of x — X'^{y, if) induces a gain in e and is stated as follows. 

Lemma 2 Let u G 5[+oo,3c?] and V{x,y,r]) a -valued function in S[+oo,3d], then 

(k; V{x, y, r])-{x- X^{y, r]))u; 8) = iel^ {k; L{k; 8; V)u; 8) (11) 

where 

[i(K;8;F)M](a;,y,7?) := div„ L(x, y, r;):y(F«(y, 77); 8)t" V(x, y, ry)! • (12) 
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It remains to remove the possible polynomial growth in 77 which follows from the relation 

where A{F'^{y, r/); 0) stands for the S[Q; 2d] matrix given by 

A{F; 6) = [C^ - lA^e] y{F; Q)~^ = [C^ - lA^O] [D^ - iB^e] . (13) 
Lemma 3 Let u E S[+oo, 3d] and V{x, y, 77) a C^-valued function in 5[+oo, 3d], then 
(k; V{x, y, 77) • 77U; 6) = -ieJ^ (k; L'{k] 6; V)u; 9) 

where 

[L'{K;Q;V)u]{x,y,'q) := divj^ [it(a:, y, 77)1/(0;, y, 77)] +div^ [u{x,y,r])A{F''{y,r]);Q)'^V{x,y,r])] . 

4 Continuity 

In this section, we will prove an i^-boundedness result for our FIOs analogous to the 
Calderon-Vaillancourt Theorem (see [S]) for pseudodifferential operators. 

Theorem 2 Let u E S[0;3d] be a symbol and k a canonical transformation of class B, then 
T^(k;w;0) can be extended in a unique way to a linear bounded operator L'^{M.'^) L'^{M.'^) 
and there exists C{k; Q) > such that 

Qj=0,l 

ft =0,1 

In the situation of a compactly supported symbol independent of x already mentioned in 
[T5] and [3] , it is rather straightforward to adapt the proof of Proposition 5 in [15] to get the 
following result. 

Theorem 3 Let u E Cj?°(M^'^; C) be a compactly supported symbol in {y, j]) and n a canonical 
transformation, thenI'^{K;u] &) can be extended in a unique way to a linear bounded operator 
L'^iR'^) L'^iR'^) and there exists C{n; 6) > such that 

\\I'{K-u-mL^^L^^C{K-Q) J2 \W\oo- (15) 

Proof of Theorem [2] We first note that the scalings of the third point of Remark [3] 
indicate that it is enough to consider the case e = 1. We will then use the strategy of Hwang 
in [H]. We have 

{^\l\K;u;e)^) = [ e'*"(-'y^''^«'u(a;,y,77)^^(7/)d7/d.Tdr; 

e'[*"(^'«'"'®)-«-«lu(a;, y, 77)VK^(J^" V)($)dCd7/da;dr;. 
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(27r)3'^/2 

A straightforward computation shows that 

(V, + A{F-{y, 77); e)V„) [<i>-{x, y, 77; 9) - ^.y] + 

where A{F;Q) is defined in (|13p . Thus, one can integrate by parts d times to write the 
integrand in (-0 | 1^{k; u; Q)(p) as 



Jl>S>'^ix,y,7j;0)~i.y] 



^ 1 + tdy^ + ^ Eti dy. ° AjkiF^jy, y); 9) 



u{x,y,ri)ip{x){T V)(C) 



^ en*''(^,...;e)-,..];,^^^(y^ ^. 0) ^gpg,^^^^ -Q (^^^^^^.^^^.)),^,^, (^-V)(0 



ft =0,1 
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with hfj^s € S'[0, 2(i]. Hence, by Fubini theorem 



l7K<i 



where we have set 



{2Try'i/^ jj-^ (1 - + + ^ 

Now, we will intcrprctc those functions as the short-time Fourier transform (or cross- Wigner 
distribution) of functions and use the results and notations of Appendix |BJ Precisely, we 
have (see Definition [5] and formula (|22p) 

^f^,{X,E;K;e)^Vg^^m-X,E) and .i.s{y,v)^j^^^Vp,[^-'^]{y,~7j) 

where 

d 

50^(x,X,S;e) = e-^^-«-(9^a»(x+X,.i-i(X,S)) and ^^(e) = J] -— — — ^. (16) 
Thus, 

{ij \l\n;u-Q)^) ^ IIWIIl~-II*/37IU=-II*5|U= 

ft =0,1 

< E IIWIlL-c'(e)||5«5,^9>|U~||v^lU.-^i^lb^ 

Qj=0,l; ft=0,l ^ ' 



^ 7r'*/2C(e) 



aj=0,l; ft-=0,l 

which closes the proof. □ 



5 From Frozen to Thawed Gaussians 

As presented before, the complex matrix does not depend on {x,y,rf). For applications, 
it might be useful to allow such a dependence. Obviously, as 9 should control the spreading 
of the Gaussian in a; around the phase-space point K(y,r]), it is rather unadequate to allow 
for ^-dependence (otherwise it is not a Gaussian anymore). Thus, we will treat the situation 
where 0(y, rj) only depends on y and 77. The following definition gives a precise idea of which 
dependence is allowed. 

Definition 6 (Admissible spreading) Let Q be a function on R"^ X with values in the 
set of complex symmetric matrices of positive definite real part. We will say that Q is an 
admissible spreading matrix (Q C) if Q € S[Q;2d] and there exists a constant positive 
definite real symmetric matrix Qq such that 

V(?/, 77) e M'' X K'', VFgR'', V ■m{y,r])V -QoV. 
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The definition of FIOs for those general admissible spreading matrices is mostly the same. 



Definition 7 (FI02) Let Q be an admissible spreading matrix in C. For u G S'[+cx3;3o?] 
and Lp e 5(M^; C), we define the action on Lp of the semiclassical FIO associated with k 
with symbol u as the oscillatory integral 



[T{K.;u;Q)ip]{x) 



1 



(27re)^ 



^i-f^ix, y,nmy,v)) [u{x, y, ri)ip{y)] dr]dy. 



Let us now justify that Theorems [T] and [2] still holds for those FIOs. The key idea, as 
presented in |16| , is to observe that the non-constant spreading part could be driven into the 
symbol without any damage. More precisely, we will split the Gaussian into 



exp 



1 



&{y,ri):. 



exp 



1 

2e 



eo 



X exp 



1 



en 



Thus, we have 



I=(/t;M;e) =X^(K;w^[e];eo) 



with the new symbol 



v'[Q]ix,y,v) = "(a;,y,'7)exp ( -^{x - X'^{y,rj)) 



(a--X«(2/,r;)) 



which is 5[+oo,3(i] if u does (the polynomial growth appearing when differentiating are 
compensated by the Gaussian decay with spreading still eo/2). As for the extension of 
Theorem [21 it is enough to notice that the scalings of the third point of Remark [3] implies a 
stronger result where each partial derivative of ([T^ induces a gain of ^/e which compensates 
the ^/e ^ that appears when differentiating the symbol [Q] . 



6 Application to Evolution Equations 

Theorem 4 Let e~i^ * be the propagator defined by the time- dependent Schrodinger equa- 
tion ([T]) on the time-interval [— T, T] with subquadratic potential V G C°°(K'',R) and let 
Q G C^{\—T,T]);C) be a time- dependent admissible spreading matrix such that Q{t = 0) is 
constant in {y,ri). Then 



sup 

te[-T,T] 



N 



n=0 



(17) 



where k* = (X" , ) and the Un{t,y,r]) are uniquely given as 

• the Hamiltonian flow associated with h{x,£,) — ^ — h V{x) and 

• the solutions of the Cauchy-problems 



—Unit,y,r]) 



1 



Un{t,y,r]) 



with initial conditions 



^L„^k['d]ukit,y,T]) 



k=0 



uo{0,y,ri) 

Un{0,y,v) 



The matrix y function of F'^ , the Jacobian of , is defined in ^ and the Ln_k's are 
differential operators whose coefficients are functions of derivatives of the Hamiltonian 
h and of the admissible spreading matrix 0. 
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Moreover, for the Ehrenfest timescales T[e) = Ct | In e|, the constant C{T) of (jl7p becomes 
C'£~'' with p arbitrary small if Ct is taken small enough. 

Remarks 4 

1. The choice the action associated to is made by continuity in time with the initial 
condition S'^ = and turns to be exactly the classical action integral defined by 



ds 



X-\y,Tj)-E^\y,Tj)-h{K%y,r^)) 



ds. 



2. The in,fc[Q] 's can be obtained recursively by the process explained in the proof. More- 
over, in the case where & does not depend on {y,r/), L„,fc[8] only depends on n — k 
and involves only derivatives of h of order between n — k and 2{n — k), their formulae 
can be made rather explicit using the annihilation/creation point of view as presented 
in [17] and fl6l. 



Proof When O does not depend on (y, 77), the way to obtain recursive equations for the 
u„'s was already explained in detail in [13j and [5] and strongly relies on Lemma [2] Indeed, 
one easily gets, with u'^ 



e-u 



is—T{n';u''-Q{t)) = I" 
H^I^{K*-u^-e{t)) 



K ; le—u 
dt 



dt 



2 



V{x) 



u^]Q{t) 



As if^''' 

at 



\J x'^^ P and Aj;$'' are polynomial in x, one can remove this dependence using 
Lemma [2] and the identity x = X'^* + {x — X'^'). To treat the V{x) term, one uses Taylor 
formula at point X*^* to the (2N + Vf^ order and treat the polynomial appearing as before. 
As for the remainder, iterative applications of Lemma O show that the gain in e is at least 
with a symbol still dependent on x. Thus, we obtain 



ie— - 
dt 



N+l 



I'{k';u'- Bit)) = ©W) + s^+^TiK';v'j^_,^; 9(0) (18) 



where the symbols on the right hand side are explicit functions of the u„'s. More explicitly, 
we have 



dt dt 



ho K 



uo 



dt 



(0 ^ n sC iV) 



fe=0 



where the L„_fc[0]'s are linear differential operators with coefficients in S[0,2d] (more pre- 
cisely, they are polynomials in derivatives of the flow k* and in derivatives, of at least second 
order, of the potential V{x)). The cancellation of the prefactor in vq (we do not want uq to 
vanish identically) shows that k* must be the flow of h. Next, cancellations of the w„'s for 
1 ^ 71 ^ + 1 give the transport equations for the w„'s whose solutions are easily seen to 
be bounded. The final step is then to check that the remaining symbol 'y|/_|_2 is bounded so 
that Theorem [5] applies. 

For non-constant admissible spreading matrix, the proof follows the same lines. However, 
we need to reformulate Lemma [2] because we have now 

{v ■ V„) <i>^ = v y{F^- e)(.T - X'^) + '-{x - A«) • [{v ■ V,) 6] {x - A«) 



so that the lemma involved is now 
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Lemma 4 Let u G S'[+oo, 3(i] and V{x, y, if) a 'C'^ -valued function in 5[+oo, Sd], then 
I' {k- V{x, y, 77) • {x - X^iy, r,))u; 9) = (19) 
iel' {k; L{k; 9; V)u; 9) + ^I" (k; {x - X") ■ yiF""; 9)^"^ [{V ■ V,,) 9] {x - X^)u; 9) 

where L{k;Q;V) is defined in (jl2p . 

Thus, one removes one (x — X'^) (gaining one order in e) to the price of adding a quadratic 
term in (x— X"). To summarize, either we gain in orders of £ or either we produce polynomials 
in {x — X'^) of higher degree. Thus, a straightforward induction shows that there exist linear 
differential operators in,Q[9] and Lk.a,p[Q] such that 

N+l 

I^n-Vo,{x,y,'n)(x-X^{y,7^)Yu-Q) = ^ (k; L„,„[9](147/); 9) 

+ e"'!' {K-{x-X'^fLk.o.,pm{V^u)-Q) 

fc+M J.Ar+2 

To estimate the remainder, one combines the polynomial [x — X^)^ with the Gaussian 
part of exp to obtain the appropriate gain in orders of e (see Section [5] above and 

Section 7 of [fe]): 

Lemma 5 Let u £ S'[+oo,3d], then 

\\T {K-{x-X-fu-Q)\\^Cpme\^^'^ J2 \\d^d^d:'^u\U 

Qj=0,l 

Finally, for the Ehrenfest timescales, it is enough to see, using estimates from [I] or [3], 
that the explicit forms of the u„'s for n 2 and thus of the symbol in the remainder 
involve derivatives of the potential (of at least second order) which are assumed to be globally 
bounded and derivatives of the linearized flow F'^ which have uniform bounds C'£~f as 
already stated in Propostion 1 of [T7] . □ 

Remark 5 The fact that our class of FIOs is not stable under the action on the left of 
pseudodifjerential operators (some of them, though bounded, can not be represented as FIOs 
with K = id and a bounded symbol) makes the generalization of this Theorem to general 
Hamiltonians rather technical. However, following the lines of Section 5 in [TB] , it should not 
be out of reach. 



A Oscillatory Integral with Complex Phase 

We present here the standard machinery of oscillatory integrals. For the definition of expres- 
sions like 

^_ I et*(-'^-'')a(a;,y,77)d77dy, (20) 

which have no sense as an ordinary Lebcsgue-integral, two approaches can be taken. First, 
one can choose a function a G S{R'^) with a{0) = 1 and set 

(I2nD:= lim [ a(77/A)ei*(^^«''')a(a:,y,77)d77dj/. 

To show the independence of the function cr a second technique is required. Under suitable 
conditions on the phase fimction, sec for instance |14| , the operator 



L 



1+ |B(y,r?)V(y,^)$(x,y,77)|2 



1 - ieB{y, r7)V(.y^,,)$(a;, y, 77) • B{y, v)'^{y^r,) 
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provides decay in rj by partial integrations, i.e. 

where is the symmetric of L defined by 

{Lip)iy)^P{y)dy - / cp{y)iL^^){y)dy, V^, ^ G 5(M'^) 



Hence an alternative definition is provided by 

For the special case of the phase function the operator L reads 

l + i£?/- [Vy+^(y,7?;e)V,,] 



1 + 



where A{y, rj] 9) is defined in (fT^ . and provides the expected decay in the 77- variable. As for 
the decay in the y variable, it comes from the Schwartz class of ip in a(x, y, 77) — u(x, y, ri)ip{y). 

B The Short-Time Fourier Transform 

We begin with a definition. 

Definition 8 For f and g two functions in the Schwartz space S{M.'^), we define their short- 
time Fourier transform as the function on x 



VgifKy^v) = J^^ e-^^-^f{x)g{x - y)dx. 

Lemma 6 The bilinear operator {f,g) E 5(]R'') x iS(]R'^') 1-^ Vg[f] extends by continuity to an 
operator L^(R'') x L'^{R'^) L^iW^ x W^) and 

\\vM\l-^ = II.9IU-II./IIl- 

Proof We use Parseval formula in rj and Fubini theorem to get 

llK,[/]||i. = ||^-V,[/]||i. = = / \fix)\'( [ \g{x^y)\'dy)dx. 

Hence, the result. □ 

The reader interested in more properties of this transform should consult [8] . Anyway, we 
extend here the definition of short-time Fourier transform for window functions g Gaussian 
in X but also dependent on the variables y and 77, more precisely of the form 

g{x, y, ?/) = G{x)g{x, y, 77) (21) 

with G{x) = exp(— a; • Qx/2) and g e 5[(+oo,0); {d,2d)], by the formula 

1 

(270"^/2 



^9 if] (2/: V) = / - y, y, rj)dx. (22) 



The preceeding lemma also extends to this situation and states the following. 
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Lemma 7 If g is of the form (|2ip . then the linear operator f >—>■ Vg[f] extends by continuity 
to L2(Mrf) ^ L2^R'i X M'') and there exists a constant C[e] > such that, for any f € L'^{M.'^), 



v,[f] ^ cm I w^rgh^ 1 ii/iu^ 



I a, =0,l 



where the sum runs on all multi-indices a = (ai, . . . , ad) such that any aj ~ or I. 



Proof We write 



Vg[f] as the integral in (.t^, x^, y, rf) of 



L2 



Integrating by parts with respect to and x"^ turns the integrand into sum of terms 



-{Tf){e)g2ix'-y,y,v) 



n(l + ^(ef-^,)) 



with gi and g2 of the form (|2ip and whose integral reads 



"■(-'--'^Fp^ 7?)^^, rj)g2ix' - y, y, 77)51 (xi - y, y, r,)dx'dx^dyd7j 



where po stands for the function defined in ^TE\i . Moreover, the Gaussian parts of gi and 
52 combine into 



exp 



so that 



1/2 



x^ + x'^ 



— ■ — h i{x — X ) ■ ^^yB 



Vg[f] is a finite sum of integral in {x^, x'^, y, rj) of terms of the form 



L2 



Mx\v)Mx^v)y''e-\<-^^^'''-y\\x' -xye~'^\(^^'^^^^ 

with Vi,V'2 e L2(R2'^;C), f/'s e 5'[(+oo,0), (2d,2d)] and a,/3 € N'^. As estimating those 
integral, up to raising a and (3 to take into account the polynomial growth of ips, we end up 
with, after the change of variables {x^,x'^) = (x,x — z). 



II^sIIl- I / \y\e 

^ \\H\l 



a, -|(!Re)i/ = a|^ 



dy 



\y |e 



m^y''y\'dy] f / |z''|e-^i(^Q)^^^^i^dz)|K,|U.||v.2|| 



We conclude using Ij^AilU^ = II^2||l2 = 7r'^/^||/||i2 (because of UpolU^ = n'^/^). □ 
Acknowledgements : The author would like to thank C. Fermanian Kammerer for fruitful 
discussions and J. Le Rousseau for indicating him the strategy of Hwang. 
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